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I. INTRODUCTION 

One of the main aims of quantum metrology is to find 
the ultimate precision bound on the estimation of a phys¬ 
ical parameter encoded in a quantum state. Of special 
interest are parameters that cannot be measured directly, 
since they do not correspond to observables of the system. 
However, they can be estimated by finding an appropri¬ 
ate measurement strategy. The estimation also involves 
choosing an estimator e which maps the set of the mea¬ 
surement results onto the set of possible parameters. The 
ultimate precision limit is given by the quantum Cramer- 
Rao bound mm which gives a lower bound on the mean 
squared error of any locally unbiased estimator e. The lo¬ 
cal unbiasedness means that in the limit where the num¬ 
ber of measurements goes to infinity, the value of the esti¬ 
mator converges to the real value of the parameter. The 
bound is given by the number of measurements taken on 
the identical copies of the state p{e) and a quantity H{e) 
called the quantum Fisher information. Higher precision 
is achieved by increasing the number of measurements 
and maximizing the quantum Fisher information. Cal¬ 
culating the quantum Fisher information thus gives us 
an idea of how well we can estimate the parameter when 
only a fixed amount of measurements are available. This 
technique has been applied, for example, in large inter¬ 
ferometers like VIRGO [3] and LIGO [3] assigned to mea¬ 
sure gravitational waves, or a current proposal [5] of mea¬ 
suring gravitation waves using phonons in Bose-Einstein 
condensates, magnetometers HUH, and gravimeters 0. 

Calculating the quantum Fisher information is not al¬ 
ways an easy task. Although a general formula for the 
quantum Fisher information exists, it is written in a 
terms of the density matrix [S|. On the other hand, many 
applications use a special kind of continuous-variable sys¬ 
tems called Gaussian states, for which the description us¬ 
ing density matrices seems particularly ineffective. Gaus¬ 
sian states can be conveniently described in terms of the 
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first and the second moments of the so-called quadrature 
operators. This description is usually called the phase- 
space or the covariance matrix formalism [Sj. 

Despite of the importance and practical usage of the 
quantum Fisher information, the theory for estimation 
using Gaussian states in the phase-space formalism is far 
from complete, and only partial results are known. The 
first leap in deriving general formulae has been taken by 
Pinel et al. [TOj, who found a formula for pure states, i.e., 
for the states which are pure at point e and remain pure 
even if the e slightly changes. The same year Marian and 
Marian found the formula for the fidelity between one¬ 
mode and two-mode Gaussian states m , which allowed 
for the derivation of the general formula for the one-mode 
state |T5]. Also, Spedalieri et al. found a formula for 
the fidelity between one pure and one mixed Gaussian 
state m , from which one can derive a slightly more gen¬ 
eral formula for pure states, i.e., for the states which are 
pure at the point e but a small change in e introduces 
impurity. A different path has been followed by Mon- 
ras |14j . who connected the quantum Fisher information 
to the solution of the so-called Stein-equation. Using this 
approach, he derived the quantum Fisher information for 
a generalization of pure states called iso-thermal states, 
and a general formula for any multi-mode Gaussian state 
in terms of an infinite series. Using the previous result, 
Jiang derived a formula m for Gaussian states in ex¬ 
ponential form and simplified a known formula for pure 
states. Quite recently, Gao and Lee derived an exact 
formula m for the quantum Fisher information for the 
multi-mode Gaussian states in terms of the inverse of cer¬ 
tain tensor products, elegantly generalizing the previous 
results, however with some possible drawbacks, especially 
in the necessity of inverting relatively large matrices. 

In this article, we first introduce a phase-space descrip¬ 
tion of the Gaussian states. Then we derive an exact for¬ 
mula for two-mode Gaussian states in analogy with [12] . 
We simplify the result of Monras for multi-mode Gaus¬ 
sian states and we show that the infinite sum involved 
converges as a geometrical series. However, the series 
may not be easy to evaluate, and for that reason we cal¬ 
culate the error when only a finite number of terms of the 
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sum are taken into account. On the other hand, in the 
case when the Williamson decomposition of the covari¬ 
ance matrix is known, the infinite series can be evaluated. 
This gives a general and exact formula for the quantum 
Fisher information of any multi-mode Gaussian state in 
terms of its decomposition. Finally, we use the derived 
formula to find the optimal Gaussian states for estimat¬ 
ing the squeezing parameter of a single mode squeezing 
channel, and we demonstrate that this estimation can 
be exponentially enhanced by the initial squeezing when 
followed by an appropriate rotation. In [A] we connect 
the so-called real and the complex form of the covari¬ 
ance matrix, in [5] we study the case of pure states and 
the regularization procedure which allows us to use our 
results for generally mixed states, even for states where 
some or all modes remain pure. JO and [D] contain detailed 
proofs for some results and in m] we show a full deriva¬ 
tion of our single mode squeezing channel example. In|F| 
we provide a table of frequently used notation that are 
constant throughout the paper and appear repeatedly. 


A. Phase-space formalism of Bosonic modes and 
the Williamson decomposition 


In this section we recapitulate the phase-space de¬ 
scription of a Bosonic system, which will be particu¬ 
larly useful for the continuous parameter states known 
as Gaussian states. First, let us consider a Bosonic sys¬ 
tem with the set of annihilation and creation operators 
{a n ,al l }. We collect them into a vector of operators 
A := (ai, a 2 ,..., o|, a\, ■ ■ ■ ) T ■ Now we can write the com¬ 
mutation relation between these operators in a compact 
form 


[Ai, A]] = Kijid 


K = 


I O' 

0 -I ’ 


( 1 ) 


where id denotes the identity element of an algebra and 
I is the identity matrix. Note that AT -1 = A'f = K and 
K 2 = I. 

One way of representing a state in quantum mechan¬ 
ics is by using a density matrix p , however, for Bosonic 
systems an alternative and completely equivalent descrip¬ 
tion exists, which is particularly useful in a description of 
Gaussian states. Given a state p we define the symmetric 
characteristic function as 

X(0 = tr[p £)(£)], (2) 

where l)(£) = is the Weyl displacement operator 

with a variable of the form £ = 707 . Gaussian states 
are those whose characteristic function is, by definition, 
of Gaussian form, i.e., 

x(0 = (3) 

In analogy with classical probability theory, Gaussian 
states are completely described by the first and the sec¬ 
ond statistical moments d and ct, where the vector d and 


the positive-definite Hermitian matrix o are defined as 


tr[Mi], 

(4a) 

tr[/>{AA*, AA]}], 

(4b) 


where A A := A - d, and {-,-} denotes the anti¬ 
commutator. We call d the displacement vector and 
o the covariance matrix. Note that other authors use 
different conventions. We choose the covariance matrix 
defined by the anti-commutator of annihilation and cre¬ 
ation operators which is known as the ‘complex form’, 
while some authors define it using the correlations be¬ 
tween position and momenta operators. In our conven¬ 
tion the vacuum is represented by the identity matrix 
/, i.e., the variance of the quadrature operators oti and 
Pi are var(ii) = var(pi) = +1 (some authors define the 
vacuum variances as +1/2). This is of course only a defi¬ 
nition and does not affect any physical interpretation of 
the results. 

According to Williamson’s theorem [171419) any 
positive-definite 2 N x 2 N Hermitian matrix er can be di¬ 
agonalized using symplectic matrices, i.e. er = SDS\ 
where S is an element of a complex representation m 
of the real symplectic group Sp(2N,M), i.e., S is an ele¬ 
ment of a group isomorphic to the Sp(2N,M.), and D is 
a diagonal matrix. S and D take the form 


a 0 
0 a ’ 


A> = 


L O' 

0 L ’ 


(5) 


where S additionally satisfies the relation SKS 1 = K 
with K defined by Eq. |T]), and L = diag(Ai, ..., Ajv) is 
a diagonal matrix consisting of the so-called symplectic 
eigenvalues of a covariance matrix er. This result will be 
used through this article. 

The symplectic eigenvalues can be found by solving the 
usual eigenvalue problem for the matrix Ko. Eigenval¬ 
ues of A'er always appear in pairs. If A i is an eigenvalue 
of Ko , then also -A i is an eigenvalue of the same op¬ 
erator. The symplectic spectrum is then defined as a 
collection of the positive eigenvalues of Ko. In other 
words, A; is a symplectic eigenvalue of o if and only if 
it is positive and ±A; are the eigenvalues of the oper¬ 
ator Ko. Symplectic eigenvalues are always greater or 
equal to one and are related to the purity of the Gaus¬ 
sian state. The state is pure if and only if for all i, Aj = 1, 
and the larger the symplectic eigenvalues are, the more 
mixed the state is. Knowing this, we can say symplec¬ 
tic eigenvalues are analogous to the eigenvalues of the 
density matrix p in the density-matrix formalism. On 
the other hand, symplectic matrices S usually represent 
some form of a squeezing or an entangling operation and 
are analogous to the unitary operators in the density ma¬ 
trix formalism. Given the special form in Eq. ([5| of the 
symplectic matrices and the relation SKS^ = K , one can 
easily prove [221 that the complex form of the symplectic 
matrices forms a subgroup of the more general pseudo¬ 
unitary group U(N,N ) = {S e GL(2N,C)\SKS^ = K}. 
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For more details about the complex and the real form 
of the covariance matrix see lAl for a more detailed anal¬ 
ysis of Gaussian states see [9FT2T]. 


II. QUANTUM ESTIMATION OF TWO-MODE 
GAUSSIAN STATE 

In this section we derive an exact expression for the 
quantum Fisher information for any two-mode Gaus¬ 
sian state. There are numerous ways to compute this 
quantity, however, for the purpose of this section we 
adopt the definition via the Bures distance |22j . The 
Bures distance is a measure of distinguishability be¬ 
tween two quantum states P 12 and is defined through 

the Uhlmann fidelity [23! F(pi,p 2 ) : = (try/ \/~P\ P 2 \fPiY 
as d 2 B = 2(l - \JT(pi, p 2 ))- The quantum Fisher infor¬ 
mation which measures how well we can distinguish two 
neighboring states p e and p t+ dt is defined as a limit [24] 

H(e) = 8 lim 1 ~ V^iPe, Pe +de± 
de^O de 2 

The problem of finding the quantum Fisher information 
thus reduces to expanding the fidelity around the point e. 
As stated before, for Gaussian states the density matrix 
can be represented by a couple of the first and the second 
moments, p± = (di,ai), P 2 = (^ 2 , 02 ). In the case of 
a two-mode Gaussian state the fidelity can be written 
as [TT1 


A -ScP(tj 1 +a 2 ) 1 8d 

Hpi,P2) = - , 2 , ( 7 ) 

{Vf + VX)-^/{Vf + VX) -a 

where Sd = di-d 2 is the relative displacement and A, F, A 
denotes three determinants defined as 


A = |(7i + cr 2 |, 

(8a) 

T = \I + KaiKa 2 \, 

(8b) 

A = |(7i + K\\a 2 + K\, 

(8c) 


with K = I © -I already introduced in the previous sec¬ 
tion. We note that our definitions differs from m by the 
factor of 2, which results in the factor of 4 in Eq. ([7|. 

Let us denote an expansion of an arbitrary matrix 
around point e up to second order in de as 

M = M 0 + M\de + M 2 de 2 + 0(de 3 ), (9) 

where M = M(e + de), Mq = M(e), Mi = M(e), M 2 = 
|M(e), where dot denotes the derivative with respect 
to e. Using this notation and the definition we can 
write H(e) = -4P 2 (e), where T from now on denotes 
the fidelity between the two close states p e and p e+ de- 
The problem is that if we try to use Eqs. ^ to expand 
Eq. © directly to find the expression of P 2 , we arrive 
at a complicated expression that depends on the second 


derivatives of o. However, expressions of the quantum 
Fisher information in (2j [14] do not depend on second 
derivatives. This is true most of the time but there are 
places where the second derivatives appear. This fact is 
usually ignored in the literature and we will explain those 
scenarios later. First, to find an expression that depends 
only on first derivatives we use Williamson’s theorem o - 
SDS I to rewrite Eqs. ([s]) as 


A = \D 0 + PDP^\, 

(10a) 

r = |/ + DqPdp^\, 

(10b) 

A = |A) + K\\D + K\, 

(10c) 


where the matrix P is defined as P = S 0 1 S. Following 
Eq. P has expansion P - I + P\de + P 2 de 2 + 0(de 3 ), 
where Pi = S^Si = S'(e) _1 S'(e) and P 2 = ^q 1 ^- A useful 
property of P is that it is symplectic, i.e., PKP' = K , 
giving us conditions on the first and second derivatives, 

P] = -KPiK, (11a) 

P 2 f = -KP 2 K + KP 2 K. (lib) 

Using the expansion of the determinant, 


|Mo + M\d€ + .M^e^ + O(de^)| — |A/o|+ (trAfo ^dt 

+ -htvMo 1 M 2 + (tr Mq 1 Mi) 2 - tr [ (M 0 ~ 1 M x ) 12 ]) de 2 ) + O ( de 3 ), 

( 12 ) 

which holds for an invertible matrix M 0 , equations (JTTj) , 
and the cyclic property of the trace we remove the de¬ 
pendence on second derivatives P 2 and obtain 


A = |A)|(l6 + 8tr[Do 1 .Di]de + UtT[P 2 ] 

- ^[D^KPiDqKPi] + 8tr[P q 1 ^] (13a) 

+ 2(tr[Pp 1 P 1 ]) 2 - 2tv[(Do 1 Di) 2 ])de 2 ) + 0(de 3 ), 

T = |G|(l + tr[C'~ 1 P 0 T ) i]* + (tr[G _1 P 1 2 ] 

- tr[(C'“ 1 Pi) 2 ] - tr[(C'- 1 P 0 ^Pi) 2 ] (13b) 

+ tiiC~ 1 D 0 D 2 ] + ^(triC^DoDi]) 2 

-hr[(C- 1 P 0 7?i) 2 ])de 2 ) + O(de 3 ), 

A = |P| 2 (l + tr[P _1 Pi]de + (tr [E~ 1 D 2 ] (13c) 

+ ^(trtP- 1 ^]) 2 - UviiE-'Dif^de 2 ) + 0(de 3 ), 


where we have denoted C = I + Dq and E = Dq + K. In 
the above we assumed that there exist a Taylor expansion 
around point e in this particular form, which is not true 
for A in the case where at least one of the symplectic 
eigenvalues of a is equal to one. We address this subtle 
issue in[B] Assuming all symplectic eigenvalues are larger 
than one, we insert expressions (13) into the Uhlmann 
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fidelity Q. We derive that the zeroth order sums to 1, 
the first order vanishes, and the second order provides 
the quantum Fisher information 


H(e) 


A 2 - 8((n/T) 2 + (x/A) 2 ) 

2(^+7^-4) 


+ 2d\a Q ^ d \. 


(14) 


The denominator in the above expression actually has 
a compact expression in terms of either a) the deter¬ 
minant of the state or b) the symplectic eigenvalues, 
+ n/AO - 4 = 2(|tr| - 1) = 2(|Z?o| - 1), which helps 
with computations considerably. To derive an expression 
which depends only on the first derivatives, we insert 
Eqs. (13) into the above formula. The terms propor¬ 
tional to Z? 2 vanish giving us an alternative expression 


for the quantum Fisher information, 

H (e) = j^j- T (|i?o|(tr[P 1 2 ] -tvlD^KP.DoKP,]) 

+ \^(tr[(C- 1 Pi) 2 ]+tr[(C- 1 P 0 ^Pi) 2 ]-tr[(7- 1 P 1 2 ])j 
+ Jtr[(A) + Ky'D^Dl] + 2d\ao 1 d 1 , 

(15) 

where Pi = S^Si = S(e)~ 1 S(e), Dq = D(e), D\ = D(e), 
C = I + Dq , and di = d(e). The above formula is useful 
when we know the initial Williamson decomposition. For 
example when we are trying to estimate squeezing in the 
case of two-mode squeezed thermal state. However, in 
general finding the Williamson decomposition is not an 
easy task. That is why we find an alternative expression 
only in terms of the covariance matrix a, displacement 
d and symplectic eigenvalues. For convenience we use 
the dot- notation, where dot denotes the derivative with 
respect to e. For convenience we also denote A := Ka(e). 
The quantum Fisher information for a two-mode Gaus¬ 
sian state is given by 


P(e) = 


2(14-1) 


|H|tr[(H l A) j + y/j7D42|tr[((/+ H 2 ) 1 A) ] + 4(A 2 - A 2 )| - —^ j j 


+ 2d'a~ 1 d, (16) 


where the symplectic eigenvalues of a can be calculated 
as 


Ai, 2 = ^[^iV^^ 2 ]) 2 -^!, (17) 


and | • | denotes the determinant. Using the Williamson 
theorem one can prove that the above formula reduces 
to Eq. (15). Also, strictly speaking, the above formula is 
defined only for covariance matrices with both symplectic 
eigenvalues larger than one. However, we can use a regu¬ 
larization procedure which allows us to use this formula 
in any case. This consists of multiplying the original co- 
variance matrix a by a mixedness parameter v > 1, using 
the formula (16) to calculate the quantum Fisher infor¬ 
mation for the state va , performing the limit v —* 1, and 
adding the second derivative \(e) for every symplectic 
eigenvalue A,(e) which is equal to one. We need to add 
these second derivatives because by performing the limit 

V 2 - 

F-»lwe set the problematic terms of Eq. (161 to 


zero. However, in cases when the symplectic eigenvalue 
A i is equal to one, such terms have a non-zero contri¬ 
bution which needs to be accounted for. Altogether, we 
have 


P(e) = limP(j/cr(e)) + y Aj(e). (18) 

^ i:Ai(e)=1 


More details about the procedure can be found in[B| 

Note that in the case where the symplectic eigenvalues 
do not change with a small variation in e, i.e. Ai = A 2 = 0, 
or Ai = 0 for Ai = 1, the term depending on the symplec¬ 
tic eigenvalues vanishes. This includes the case where 
purity does not change or where the parameter of inter¬ 
est e was encoded into the initial state by a symplectic 
transformation. Although all the computations were per¬ 
formed in the complex representation of the covariance 
matrix, we can easily transform the result to the real 
representation. This is done by substituting A = iflan, 
d o~ 1 d —» where ifl is the real form equiva¬ 

lent to matrix K, ctr the real form covariance matrix and 
dp the real form displacement. For more details see [A] 

III. MULTI-MODE PARAMETER ESTIMATION 

In the previous section we derived an exact expres¬ 
sion for two-mode Gaussian states, however, in recent 
work na a general formula for the quantum Fisher in¬ 
formation was derived as a limit of a particular infinite 
series. Here we simplify the expression for the infinite 
series. We also find a bound on the remainder of the se¬ 
ries when we sum only a finite number of terms. We 
then go on to simplify an already known formula for 
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pure states and derive an exact expression for the multi- 
mode quantum Fisher information for the cases when the 
Williamson decomposition is known. 

According to |14j . the quantum Fisher information for 
a general multi-mode case which has all symplectic eigen¬ 
values larger than one can be calculated as 

H(e) = -tr[<jY] + 2d} a~ 1 d, (19) 

where Y is a solution to the so-called Stein equation [25] 

& = aYa - KYK. (20) 

If all symplectic eigenvalues of a are larger than one, the 
solution is unique and can be written as an infinite series 


Y = -^(Ka)- n (a-i)(aKy 

n =0 


( 21 ) 


Inserting (21) into (19), using (er' 1 ) = -CT _1 d'cr“ 1 , and the 
properties of matrix K we find an elegant expression for 
the multi-mode quantum Fisher information, 


1 M + 

H(e) = o E tr[( A~ n A) 2 ] + R m + 2 d'a^d, (22) 

where we use the notation A := Ka or A := ifla^ for 
the real form of the covariance matrix and Rm is the 
remainder of the series M —* oo. As we prove in 0 the 

remainder is bounded \Rm\ < with A min := 

min.j{Ai} being the smallest symplectic eigenvalue, i.e., 
the smallest positive eigenvalue of matrix A. This means 

that for A m jn > 1 we have limM_ >oc Rm = 0, the series 

converges, and we can write 

H(e) = l £ ti'[(A~ n A) 2 ] + 2dV 1 <L (23) 

2 n=l 


To calculate the quantum Fisher information for the 
states which have some eigenvalues equal to one we can, 
once again, use the regularization procedure (18). A 
small example, which has been already shown in m , is to 
consider the class of iso-thermal states given by A 2 = z/ 2 /, 
v > 1, which is equivalent to D(e) = vl. For such states 
we can easily evaluate the infinite sum (23) and derive 


H(e) = 


2(1+ z/ 2 ) 


tr^yT 1 !) 2 ] + 2 dV 1 d, 


(24) 


which for v - 1 gives a formula for pure states also de¬ 
rived in [10] . As noted in [15] , this formula can be fur¬ 
ther simplified. Differentiating A 2 = v 2 and multiplying 
each side by A~ x we obtain an anti-commutation relation 
Ar 1 A = -AAt 1 which together with A -2 = v~ 2 I gives 

H{e) = ~ 2(T7^) tr [ i2 ] + 2c?t<7 ‘ lc? ’ (25) 


meaning that for iso-thermal states with fixed displace¬ 
ment we do not need to invert the covariance matrix any¬ 
more. 


Using tr[(A n A) 2 ] 
2tr[Do n+2 A'~ 


= 2tr[(Dg n+1 /\ -n+1 Pi) 2 ] - 

l P 1 ] + tr[Do 2n Df] we can 
rewrite formula (23) for the quantum Fisher information 


l I\ DA'R- 


in terms of the elements of the symplectic matrices and 
eigenvalues, 


H(e) = tr[P 2 ] - £ tr[(D 2 - /)(D 0 "A'"P 1 ) 2 ] 


n= 1 


+ 2^ r [(^° + K) 1 ^ ) o 1 ^i] + ‘2d[<j 0 1 di, 


(26) 


where we have used the notation from Eq. (15). We can 


see that the diagonal part coincides with the two-mode 
case indicating the validity of the general multi-mode for¬ 
mula. However, in contrast to Eq. (23) the infinite series 


here converges also when symplectic eigenvalues are equal 
to one. Nevertheless, this does not mean that it is valid 
to use this formula for such states. To be more specific, 
plainly inserting Dq = I does not give the correct for¬ 
mula for pure states — this still needs to be obtained by 
the regularization procedure. The reason why it does not 
give the correct result is that in general limits limM^oo 
and lim„^i do not commute. On the other hand, the 
infinite sum in this formula leads to a geometric series 
that can be evaluated, as shown in[D] This allows us to 
derive a much more elegant, entirely general and exact 
formula for the quantum Fisher information in terms of 
elements of the Williamson decomposition a = SDS t for 
multi-mode Gausian states. Using the definitions ([5]) and 
(11aI, we find that matrix Pi has an elegant structure 


A(e) = 


R Q 
Q R 


(27) 


where R - a^a - /?t/3 is a skew-Hermitian and Q = 
cA $ — /?! a a (complex) symmetric matrix. Matrix Pi is 
actually an element of the Lie algebra associated with the 
complex form of the real symplectic group. If the diago¬ 
nalizing symplectic matrix forms a one-parameter group 
S = e Xe , where X is an element of the algebra indepen¬ 
dent of e, then Pi = A'. Inserting Eq. (27) into Eq. (26) 
and evaluating the infinite sum we derive the quantum 
Fisher information for the A-mode Gaussian state 


H(f , f (Ai-Aj ) 2 ■ 

R\ e ) 2.J \ \ i 

i,j=l AiAj 1 


12 ■ (A, + Ai) I a 


XiXj + 1 


Xi 


N 

Tj a 2 1 

i=l A i 1 


(28) 


+ 2d 1 a 1 d , 


where A^ are the symplectic eigenvalues. Strictly speak¬ 
ing, the formula is not defined for the symplectic eigen¬ 
values equal to one. Under the assumption of differen¬ 
tiability of the covariance matrix, for A,(e) = A j(e) = 1 

both (e) and y’jty(e) are expressions of type 2. 

Nevertheless, we can define these problematic points in a 
way which makes the quantum Fisher information a con¬ 
tinuous function, which also corresponds to the knowl¬ 
edge gained from a study of the states close to the 
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points of purity described in PBl For A,(e) = 1 we de- 
a 2 

fine ^y(e) := A;(e), and for Ai(e) = Aj(e) = 1 we define 

(Ai-A ,) 2 / ■, q 
AiAj-l A 6 '' u ' 

Assuming all symplectic eigenvalues are larger than 
one, we can define the Hermitian matrix i?. 


-j==j!=Rij and the symmetric matrix Q t j := -j==^= Q t j 

which allows us to write the quantum Fisher information 
in a compact form, 

H(e) =tr[M t + OQ t ]+tr[(L 2 -/)' 1 i 2 ] + 2d t cr~ 1 d. (29) 

Now we can easily interpret each term. The first part cor¬ 
responds to the change in relative orientation and squeez¬ 
ing, the second to the change in purity and the third to 
the change in displacement. Added together, they all 
contribute to the quantum Fisher information and in¬ 
crease the precision in the estimation of a parameter e. 
Note that formulas (28), (29) respectively, are actually 


+ At 


Hj 


multi-mode generalizations of Eq. (15). The form (15) 


takes when we rewrite it in terms of submatrices given 
by Eq. |27| ) is exactly the same as (28) for N - 2. The 
same holds for the one-mode formula derived in [12] and 
AT = 1, partially validating our general result. For the 
derivation of the formula see [D] for the details why we 
choose that particular definition of problematic points 
see [Bj and for the quantum Fisher information in the 
real form formalism see m 


IV. EXAMPLE 


To illustrate the use of the general formula (281 we de¬ 


rive the precision with which we can estimate the squeez¬ 
ing parameter e of a one-mode squeezing channel. Sim¬ 
ilar problems have already been studied. For example 
in [26j the vacuum and coherent probe states were con¬ 
sidered, in m it was the coherent states and the dis¬ 
placed squeezed vacuum states, which were squeezed in 
either position or momenta direction. Here we are gen¬ 
eralizing the precision bounds in those articles. We also 
note that a similar work has been done in E8] using the 
density matrix formalism for a different type of a one¬ 
mode squeezing channel S r = exp(-i|(a 2 + a^ 2 )). Here, 
however, we use an entirely general single-mode Gaussian 
state as an input state. 

We choose a one-mode squeezed rotated displaced ther¬ 
mal state as an input state, which is the most general 
single-mode state [21] . with the initial squeezing param¬ 
eter r, angle of rotation 9, and the initial displacement 
do = | do | e 1 ^. We feed this state into a one-mode squeez¬ 
ing channel, which encodes the unknown parameter e we 
are trying to estimate, leaving us with the final state, 

p t = S e D Jo R e S rPth SlR\D\ o Si , (30) 

where S r = exp(|(a 2 - a^ 2 )) is the one-mode squeezing 
operator, Rg = exp (-iOaJa) the rotation operator, and 



FIG. 1. Visualization of the distinguishability of the co- 
variance matrices with different initial rotations. The initial 
squeezing was set to r = 0.8, the initial displacement do = 0, 
and the final squeezing e = 0 (blue with full line) or e = 0.1 (or¬ 
ange with dashed line). The initial rotation from left to right 
9 = 0,f,f,^,f. Covariance matrices with 0 = ^ can be 
easily distinguished allowing for a better estimation of the 
parameter e. 


D rl = exp(dod^ _ dod ) the displacement operator. As we 
show in[Ej the quantum Fisher information is 


4A 2 

H(e)= 1 ( coslrr + sinlrr - 2cos(4d) cosh 2 rsinh 2 r) 

At + 1 


4 Idol 


+-( coslr(2r) + cos(2 (9 -</>)) sinh(2r)). 

Ai 


(31) 


Now, let us have a few notes on the derived formula. The 
result does not depend on the parameter we want to es¬ 
timate, which must be true for every encoding channel 
which forms a one-parameter unitary group 17(e) = e~ lKe , 
where K is a Hermitian operator. Also, as pointed out 
in [22], because the symplectic eigenvalue Ai > 1 is pro¬ 
portional to temperature, we can immediately see that 
although the thermality slightly enhances the estimation 
through squeezing (with the maximal enhancement by 
the factor of 2), it reduces the estimation from the dis¬ 
placement. More importantly, however, by choosing an 
appropriate rotation of an input state we can significantly 
increase the estimation precision. Without loss of gen¬ 
erality we assume r > 0. The maximal value of (31) is 


achieved when we choose to rotate for example by the 
value 9 = j , and displace in the direction ^ = f, 


4^2 4 ^ 

#max(<0 = 2 — cosh 2 (2r) + — |d' 0 | e 

A-. + 1 A\ 


2 r 


(32) 


which signifies an exponential increase in the precision of 
estimation the unknown parameter e when considering 
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an initially squeezed and rotated state. To demonstrate, 
we derive a formula for a squeezing needed to enhance 
by k orders of magnitude when assuming zero initial dis¬ 
placement, r = arcsinh^ 10 ^" 1 , which for larger k (or 
r > 1) behaves as r & 0.35 + 0.58/c. The current state- 
of-the-art [3D] achieves the squeezing around r = 1.46, 
which could hypothetically account for an improvement 
by a factor of 80 as compared to the case where the ini¬ 
tial squeezing is zero. As we show in figure [T] the reason 
why the amount of distinguishability of the two close 
states rises is because the final squeezing forces the co- 
variance matrix to turn. Also, note that although the 
initial squeezing leads to an exponential increase, the 
initial displacement contributes only quadratically. To 
conclude, the optimal Gaussian state for the estimation 
of the parameter of a one-mode squeezing channel is a 
thermal state infinitely squeezed in the angle of 6 - j 
from the direction of the squeezing channel we want 
to estimate, and infinitely displaced in the direction in 
which the squeezed state is stretched, (j) = f . For a fixed 
amount of squeezing and displacement, the optimal tem¬ 
perature which maximizes the quantum Fisher informa¬ 
tion is given by a solution of (A 2 A + \ )2 = 2 ^oL\ 2 r) • Also ’ 

because Ai = 1 + 2 nth, where n t h = tr[p t /ja'a] denotes 
the mean number of thermal bosons, and the mean total 
number of bosons in displaced squeezed states is given 

by 

n = n do + nth + (1 + 2n t / l ) sinh 2 A (33) 


where we have denoted nj := |do| the mean number 
of bosons coming from displacement, we can rewrite the 
maximal quantum Fisher information (32) as 


H„ 


2(1 + 2n - 2n r ) 2 

_( e ) = __ 

1 + 271^(1 + 71 ^) 


4n d 0 ( 1 + 2n - 


2n d 0 + Z^n-n^-nthy/T+n-n. 


do 


+n th \ 


(1 + 2 n th y 


(34) 


For a hxed total number of bosons this function achieves 
the maximum at n t h - n j 0 = 0. This means that for 
a Hnite amount of available energy, the best protocol is 
to invest it all into squeezing, which is the same result 
as in [28; [29] • The optimal probe state is then the j- 
rotated squeezed vacuum state with the quantum Fisher 
information H(e) = 2(1 + 2 n) 2 , which clearly indicates 
Heisenberg scaling. 


of positive-definite matrices and the properties of the real 
symplectic group. Although the formula is not directly 
applicable for states with pure modes, we introduced a 
regularization procedure which allows us to overcome this 
problem. Then, using a different approach, we simplified 
an already known formula for the multi-mode quantum 
Fisher information in terms of an infinite series. We also 
estimated the remainder of the series, allowing for an ef¬ 
fective numerical calculation. Using the previous result, 
we showed that for the cases when the Williamson decom¬ 
position of the covariance matrix is known, the quantum 
Fisher information for multi-mode Gaussian state can be 
computed exactly. The general multi-mode formula is 
equivalent to the known results for one-mode Gaussian 
state when setting TV = 1 and to the previously mentioned 
two-mode Gaussian states when TV = 2. However, we note 
that using the requirement of the continuity of the quan¬ 
tum Fisher information and studying the case of the pure 
states gave us a different definition for the problematic 
points than is mentioned in m- Finally we applied our 
newly gained formula to study the case of the estimation 
a squeezing parameter of a one-mode squeezing channel. 
We showed that a strategy of squeezing and rotating the 
input state can significantly improve the precision in es¬ 
timation. 

We believe the main achievement of this article is in 
the usefulness of the derived formulae. It allows for the 
study of the optimal input states of Gaussian nature EH. 
it helps predict the ultimate sensitivity of a physical de¬ 
tector’s particular implementation mm, or to analyse 
the effects of temperature on the current gravitational 
wave detector proposal [5] • It gives a limit in the estima¬ 
tion of time [32] or temperature [33]. Also, since certain 
objects called the Bogoliubov transformations are iso¬ 
morphic to the symplectic transformations, the natural 
application lies wherever these transformations appear. 
This is for example for quantum field theory in curved 
spacetime [SUSS] but also Bose-Einstein condensates [3S] 
or scattering problems m- 
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Appendix A: The complex and the real form of the 
covariance matrix 


V. CONCLUSION 

We have derived an exact formula for the quantum 
Fisher information of an arbitrary two-mode Gaussian 
state. This has been done using the definition of the in¬ 
finitesimal Bures distance, the Williamson decomposition 


In this section we describe the structure of covariance 
matrices and displacement. We introduce the real form 
covariance matrix formalism and find how covariance ma¬ 
trices in the real form transform into its complex form. 

The complex form of the covariance matrix and dis¬ 
placement was defined by equations 0 and §. From 
the definition we can observe the following structure of 











the first and second moments: 


d 


X Y' 

d 

, o = 

Y X 


and <7' = a, i.e., X' - X and Y T = Y. Using the 
Williamson theorem [IMS] we can write a = SDS 
where SI\S t = I\. Although the symplectic matrices 
S are not necessarily Hermitian, they follow the same 
structure as a which is expressed by Eq. (§. 

Construction of the real form of the covariance ma¬ 
trix is analogous to the complex form described in 
the introduction. It is usually defined with respect to 
the collection of quadrature operators Q := x © p = 
{x\, x 2 , ■ ■ ■ ,pi,f> 2 , • • The canonical commutation re¬ 
lations of these operators can be conveniently expressed 
as 


id 


U = 


0 I 
-I 0 


(A2) 


Consequently, a and p needed for R and Q from Eq. (28) 
can be expressed in the real form symplectic matrix ele¬ 
ments as 


u= -(ur + Sr+i^r-iPr), (A8a) 

P=^(ctR-S R +i^f R + i/3 R ). (A8b) 


Since all important matrices are related via this uni¬ 
tary transformation and traces and determinants are in¬ 
variant under such transformations, it is clear that every 
formula we derived can be easily rewritten in the real 
form formalism by the formal substitution a -*■ <j r and 
K -*■ ifl. On the other hand, the complex form provides a 
much more elegant structure and exposes the inner sym¬ 
metries of symplectic and covariance matrices in more 
detail. Also, it is much easier to work with K since it 
is diagonal, unitary, and Hermitian in contrast to non¬ 
diagonal skew-Hermitian matrix fl, providing much more 
convenient language. 


Properties of S2 are -H 2 = / and fl T = -f2, in contrast 
to the complex form version K. In the real form, the 
definitions of the first and second moments are 


d R = tr[/5Q] 


X 

p 


or = tr [p{Q:Q}] 


Xr Y r 
Yr Z r 


(A3a) 

(A3b) 


The real covariance matrix is symmetric, i.e. X R - Xp 
and Z R = Zp. The corresponding real symplectic ma¬ 
trices are given by a R = S R DpSp , where SflS T = fl, 
which is a defining relation of the real symplectic group 
Sp(2N,R). 

Since the change between real and complex form of 
the covariance matrix is a simple basis transformation, 
Q -*■ A, we can relate these two using the unitary matrix 

U, 


A.UQ. U-j. 


I +U 

I -il 


(A4) 


The resulting transformation between real and complex 
covariance matrices and displacement are 

d=Ud R , a=Ua R U\ (A5) 


and the transformations related to the Williamson de¬ 
composition are 

S = US R U\ D = UD r U ] = D r , K = Uinu f . (A6) 


We explicitly write the connection between real and com¬ 
plex form of symplectic matrix, 


Appendix B: Pure states and the regularization 


In this section we derive a formula for the QFI for the 
states around the points of purity, i.e., for the points 
where Ai(e) = 1, but not necessarily A i{e + de) = 1. We 
can see that formula (16) is undefined at points Ai(e) = 
1. This is the consequence that we assumed that for all 
A, A, T the Taylor expansion exist in a form given by 
Eqs. (13), which is not true when A i(e) = 1 for some 

iz{ !,¥ 

For simplicity let us study the case where e is such 
that the both Ai(e) = A 2 (e) = 1 and at the end of the 
calculation it will be clear how various cases work. For 
such states we have D - I + Dide + D 2 de 2 + 0(de 3 ). From 
Eqs. (10) we can see that A = 0 and A = T and the 
Uhlmann fidelity ([T]) simplifies to 


(.Pet Pe+de) ~ 


n/A ' 


(Bl) 


where for simplicity we have omitted the part consist¬ 
ing of displacement. This coincides with the general for¬ 
mula m for fidelity between one pure and one pure or 
mixed Gaussian state, which is exactly our case. Now, 
because Ai(e) = 1 and the symplectic eigenvalue cannot 
fall below 1, we have that either A i(e) =0 or the deriva¬ 
tive does not exist. If it exists and is, for example, pos¬ 
itive, then for de < 0 we have A(e + de) = 1 + A i(e)de + 
0(de 2 ) < 1 for small enough de, which cannot be true. 
From now on we assume that map a : e —*■ er(e) is dif¬ 
ferentiable, thus Di = 0 and D = I + D 2 de 2 + 0(de 3 ). The 
simplified formula ( |B1| ) for fidelity gives T (p e , p e +de) = 
1 - !(tr[P 2 ] - tr[(AWi) 2 ] + 2tr[£> 2 ])de 2 + 0(de 3 ), leading 
to the QFI, 


u R Pr 


Re [a + f}\ - Im [a - /3] 

.7 R $R_ 


Im [a + P] Re \a - P] 


(A7) 


H(e) = ^(tr[P 2 ] -trfCAPi) 2 ]) +tr[£> 2 ]. (B2) 






















9 


Since D 2 = \D, we have tr[Z? 2 ] = We can rewrite 

the above equation in terms of A = Ka, 

H(e) = J(2tr [fir'A] - tr^A^i) 2 ]). (B3) 

If the state remains pure with a small variation in e, 
then Xi = 0 and the above formula is equivalent to the 
pure state formula H{e) - jtr[(c7 _1 c7) 2 ], which can be 
also derived from the equation for iso-thermal states ( [24| . 
We see that although general mixed-state formulae do 
not depend on the second derivative, formulae ( |B2[ ) and 
(B3| do. But we can show that we achieve the same 


result from the mixed state formula when requiring the 
continuity of the QFI. If we assume that a : e —*■ a(e ) € 
C^ 2 \ i.e., the second derivative exists everywhere and is 
continuous, then the QFI is continuous everywhere apart 
from the points where it is undefined. To make the QFI 
a continuous function we define the problematic points 
given by Ai(e) = 1 as limits 


(e) := lim X ^ + de ^ 


A 2 - 1 


de-*o A i(e + de) 2 - 1 


- Aj(e), 


(B4) 


since A^ has expansion A i(e + de) = 1+ |A i(e)de 2 + 0(de 3 ). 
Now we see that the our definition of problematic val¬ 
ues using the second deriv ative corresponds to the actual 

value of the QFI given by (B2). On the other hand, note 

^ / \0. 

:= 0, which 


that authors of [12] choose rather A ^^ 2 _ 1 
would correspond to using an exac t for mula for states 
which remain pure (setting A = 0 in (B2)), but may lead 
to discontinuities even for the smooth functions er(e). 
Nevertheless, to decide which convention is the right to 
choose, one should rather examine the validity of Cramer- 
Rao bound itself for such cases. 

Similarly to definition (B4), for the multi-mode for¬ 


mula (281 and A;(e) = A j(e) = 1 we define 


(A,-A,) 
AjA^ ~ 1 


2 00 := lim = 0 , 


de-*- o A i(c + de)Xj(e + dc') — 1 


(B5) 

Now we describe the regularization process which al¬ 
lows us to derive the correct value of the quantum Fisher 
information. In fact, all (15l,(16),(23),(261 have prob¬ 


lems when there exists at least one symplectic eigenvalue 
equal to one. The basic idea is that instead of comput¬ 
ing H{e) H(a(e)) directly we calculate the quantum 
Fisher information for the regularized state H(vc ’’(e)), 
v > 1 and at the end of the computation we perform 
the limit v -*■ 1. This method must always give a finite 
result because by doing that we must arrive at an ex¬ 
pression similar to Eq. ( [28] ). This becomes clear when 
we look at the derivation in [D] Also, the result it gives 
is certainly correct for any point where the formula is 
defined, we just need to make sure that this method 
gives the same definitions (B4),(B5) for the problematic 


points. For Ai(e) = A j(e) = 1 and differentiable a we have 


A;(e) = Aj(e) = 0 and 


OMe )) 2 


= lim ■ 


0 


v ->i (vX i(e)) 2 - 1 v 2 - 1 


lim 


= lim ■ 


0 


(uX t (e) - uXj(e)f 
vXi(e)vXj(e) - 1 v 2 - 1 


= 0, 

(B6a) 

= 0. 

(B6b) 


Although the second equation gives the same definition 
as required by Eq. (B5), we can see that the first equation 


does not give the same definition as Eq. (B4) because the 


problematic term here is set to zero by the regularization 
procedure. As a result, for every e such that A i(e) = 1 
we need to add additional Ai(e) to the regularized ver¬ 
sion of the quantum Fisher information. Therefore, we 
have shown that we can calculate the quantum Fisher 
information at the points of purity as a limit 

H(e) = lim e))+ZA i (e). (B7) 

i:Ai(e) = l 


Appendix C: Estimation of the remainder in the 
multi-mode formula 


Here we prove the bound on the remainder of the gen¬ 
eral multi-mode formula. We consider the Williamson de¬ 


composition cr = SDS^. An element of the sum Eq. (22) 
can be written as 


= tr [A- n AA~ n A] = tr[D~ n BD- n B], 


(Cl) 


where B = 5'’I'A(5^) 1 K n 1 . We can derive the inequal¬ 
ities 


a n = 


I — 

< \n \n 
k,l A k A l 


BkiBik 


I — 

^ \n \n 
k,l A k A l 


\Bu\ \Bik\ 


1 


X! \Bki \I Bik 

k,l 


\ 2 n 
A mi 


At(ei b mI 2 ) 

A min V k,l / 


(C2) 


A 2, l 


where the last inequality is the Cauchy-Schwarz inequal¬ 
ity between and Bji considered as vectors with 4iV 2 
entries where N is number of modes, A m i n := min,; {A,)} is 
the smallest symplectic eigenvalue. Defining the Hermi- 
tian matrix C := AKS we have 

tr[(AA) 2 ] = tr [&DCD] = £ \C kl \ 2 X k X t 

k,l 

>AL n tr[C'tC] = AL n tr [B^B]. 
Combining (C2) and ( |C3[ ) gives 

|a„|<tr[(AA) 2 ]A m 2 r 2 . 

For A m i n > 1 we can estimate the remainder, 

tr[(AA) 2 ] - x _ 2n _ 2 tr[(AA) 2 ] 


(C3) 


(C4) 


\Rm\ < 


EA 

n=M +1 


-2n-2 

min 


2A 2 f n +2 (A 2 in -l) 


(C5) 
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Appendix D: Derivation of the exact formula for a 
multi-mode Gaussian state. 


Here we derive the general formula for the QFI given by 
Eq. (281. Denoting D oii := (D 0 ) u , we start by evaluating 


an infinite sum from the equation (26): 
gtr [(Dg-j)(Do"tf-"Pi) 2 ] 


I> [{D 2 0 - I^D-^KPtfYtvYpl - /)(D 0 - 2 ”P 1 ) 2 ] 

n=l 

E( p L-i)( ^ D ^D 0j3 y 2n+1 \KP l ) l3 {KP 1 ) jl 

i,j V n=l / 


yields 


N (Ai — Ay)~ 2 (Ai + Aj) 2 ,2 


= e 7 y: iyyy'; 

_ 2 A^ Aji JL A^ Aj u JL 


»,J=1 '"‘A? 

JV V 2 , 

+ E rv - + 2d'(7 _1 d. 


(D4) 


,tlA?-l 


The definitions of expressions for problematic points 
Aj (e) = Aj(e) = 1 come from the assumption of the conti¬ 
nuity of the QFI and correspond to the knowledge gained 
in IB 


Appendix E: Calculating the example 


+ (D'L - 1)( t^ouDojjr^yPihiPih 

=Et^Y 1 T ( (KPi) a (.KPda + DouDojj (Pi) a (Pi)ji). 

(Dl) 

We combine this expression with the first part tr[P 2 ] in 

p6|: 


tr[/?]-£(-) = E 


D U D w - 1) 


D 2 D 2 

n=l i,j ^Oii^Ojj 


(Pl)ij(Pl)ji 


D()liDojj (Dq,,; l) 

D L D 'ln - 1 


(D2) 


(KPMKPJji. 


Here, to illustrate how the general formula (28) works, 


we calculate the quantum Fisher information for an ex¬ 
ample of a one-mode Gaussian state. For TV = 1 the 


general formula (28) has a compact form 




ic 


Ai 


ill 


A?-l 


+ 2dV 1 d, (El) 


where Qn = a$ - 8a. 

We consider a task of estimating a squeezing parameter 
of a squeezing channel, with a squeezed rotated displaced 
thermal state as an input state. In the complex form of 
the covariance matrix formalism, the one-mode squeezing 
operator with a squeezing parameter r, and the rotation 
operator via angle 6 are 


Now we use definitions (|5|) , (|27|) , Pi = S^S = KS^KS, 
K = I © -I and D 0 = diag(Ai, ..., Xn, Ai, .. ., Ajv) to 
rewrite the expression in terms of the symplectic eigen¬ 
values and sub-matrices R and Q 1 


, n A 2 (A 2 -l)-AiA 7 (A|-l) 

tr [ P i ]-E(')=E 2 J \2 \ ‘i MRijRji) 


+2 


n= 1 i,j =1 Aj A j 1 

A 2 (A 2 -l) + A t A,(A 2 -l)^ ^ 


Y-l 


R e(QijQji) 


N Aj(Aj-Aj) D t D D n, , 0 Aj(Aj+Aj) D 


= y 2 

ijtl A,; Xj - 1 


Re (Ri jRji )+2 


XiXj + 1 


R e(QijQji)- 


(D3) 


Re(QijQji) and Re(R l:) Rji) are both symmetric. This 
is why we decompose parts consisting of symplectic 
eigenvalues into its symmetric and antisymmetric parts: 

0 A 7 (Aj-A j ) Xj o A, (Aj+Aj) (Aj + A,)~ 

- AiAj-l A,Aj-l XiXj-1’ A, A j +1 AiAj + 1 

A 2 -A 2 

A * A ^i and the antisymmetric part vanishes under the 


sum. In addition, from identity (11a) we have R? = 
-R and Q T = Q giving us R e(RijRji) = ~\Rij\ , 
R e (QijQji) = |Q»j| • The diagonal part is simply 
|tr[(Do + K) 1 Pq 1 Pi] ~ Ei^i yfci- Combining all these 


S r = 


cosh r 
- sinhr 


- sinh r 
coshr 


Pe 



0 

e+ie 


(E2) 


First, we create an input state by squeezing, rotating, 
and displacing a thermal state. We obtain a state given 
by its first and the second moments, 

do, a 0 = R e S r DS}.Rl, (E3) 


where do = (do,do) T is the initial displacement, do the 
initial covariance matrix, and D = diag(Ai, Ai) is the co- 
variance matrix of a thermal state. For a harmonic oscil¬ 
lator with frequency to we have Ai = coth(^|Q. Now we 
feed the prepared state into the channel which we con¬ 
sider to be again a simple squeezing operation with the 
unknown squeezing parameter e we want to estimate. We 
obtain the final state, 

d = S t d 0 , a = S e RgS r DStRlsl. (E4) 


It is clear from the construction that the diagonaliz¬ 
ing sy mpl ectic matrix needed to compute Qn in for¬ 
mula (El) is S = S e RgS r . Considering definitions ^ 
and (|E2|) we derive 


a = e l8 cosh e cosh r + e +ie sinh e sinh r, (E5a) 

8 = e~ l6 cosh e sinh r + e +l8 sinh e cosh r. (E5b) 
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For the part containing displacement we need to compute 
the inverse of the covariance matrix, cr -1 = KSD^S'K, 
and the change of disp lacement, d = S e d 0 . The middle 
part of expression (El) is simply zero, because Ai is in- 


Appendix F: Table of frequently used notation 


dependent of e. In other words, the purity of the system 
remains the same. Now we have everything prepared to 
calculate the quantum Fisher information, 


4A 2 

H(e)= 1 ( cosh 4 r + sinhr - 2cos(40) cosh 2 rsinh 2 r) 

Xf + 1 


4 \dn 


-(cosh(2r) + cos(2(0 - </>)) sinh(2r)). 


K 


S 


a, 13 

D,Dq,L 


Pi 

R,Q 

T t 


(E6) 


A constant matrix defining 

a complex represenation of the real symplectic 

group defined in Eq. 0 - 

The displacement vector, Eq. Q. 

The covariance matrix, Eq. Q. 

A := ATex, a multiple of the two previously 
mentioned matrices. 

The symplectic matrix from the Williamson 
decomposition of ex, Eq. 

Submatrices of the matrix S, Eq. ([5|. 

The diagonal matrices consisting of 
the symplectic eigenvalues, Eqs. 

The derivative with respect to the parameter 
we want to estimate. 

Pi := S' -1 .!?, an element of the Lie algebra 


associated with the symplectic group,(11a),(271 
Submatrices of the matrix Pi , Eq. (271. 
Complex conjugate, transpose, 
conjugate transpose respectively. 
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